Direct experimental evidence of a growing length scale 
accompanying the glass transition 
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Understanding glass formation is a challenge because the existence of a true glass 
state, distinct from liquid and solid, remains elusive: Glasses are liquids that have 
become too viscous to flow. An old idea, as yet unproven experimentally, is that 
the dynamics becomes sluggish as the glass transition approaches because increasingly 
larger regions of the material have to move simultaneously to allow flow. We introduce 
new multipoint dynamical susceptibilities to estimate quantitatively the size of these 
regions and provide direct experimental evidence that the glass formation of molecular 
liquids and colloidal suspensions is accompanied by growing dynamic correlation length 
scales. 



Why does the viscosity of glass-forming liquids increase 
so dramatically when approaching the glass transition? 
Despite decades of research a clear explanation of this 
phenomenon, common to materials as diverse as molec- 
ular glasses, polymers, and colloids is still lacking 
The conundrum is that the static structure of a glass 
is indistinguishable from that of the corresponding liq- 
uid, with no sign of increasing static correlation length 
scales accompanying the glass transition. Numerical sim- 
ulations performed well above the glass temperature, 
T g , reveal instead the existence of a growing dynamic 
length scale 0, EL IE IE associated to dynamic hetero- 
geneities @ • Experiments jl, IE 0, 0, 0] have indi- 
rectly suggested a characteristic length scale of about 5 
to 20 molecular diameters at T g , but its time and tem- 
perature dependencies, which are crucial for relating this 
finding to the glass transition, were not determined. 

We present quantitative experimental evidence that 
glass formation in molecular liquids and colloids is ac- 
companied by at least one growing dynamic length scale. 
We introduce experimentally accessible multipoint dy- 
namic susceptibilities that quantify the correlated nature 
of the dynamics in glass formers. Because these measure- 
ments can be made using a wide variety of techniques in 
vastly different materials, a detailed characterization of 
the microscopic mechanisms governing the formation of 
amorphous glassy states becomes possible. 

Supercooled liquids are believed to exhibit spatially 
heterogeneous dynamics over le ngth scales that grow 
when approaching the glass state [II [IE Ell EE]- This het- 
erogeneity implies the existence of significant fluctuations 
of the dynamics because the number of independently 
relaxing regions is reduced. Numerical simulations have 
focused on a "four-point" dynamic susceptibility X4(£), 
which quantifies the amplitude of spontaneous fluctua- 
tions around the average dynamics [E EL IE IE 0- The lat- 
ter is usually measured through ensemble-averaged cor- 
relators, F{t) = (SA(t)SA(0)) = (C(t)), where 8A(t) = 



A(t) — (A) represents the spontaneous fluctuation of an 
observable A(i), such as the density. Dynamic corre- 
lation leads to large fluctuations of C(t), measured by 
Xi(t) = N(6C 2 (t)), where N is the number of particles 
in the system. The susceptibility Xi(^) typically presents 
a nonmonotonic time dependence with a peak centered at 
the liquid's relaxation time 01 • The height of this peak 
is proportional to the volume within which correlated mo- 
tion takes place 0, [E 0, 0] . Unfortunately, numerical 
findings are limited to short timescales (~ 10~ 7 s) and 
temperatures far above T g . Experimentally, detecting 
spontaneous fluctuations of dynamic correlators remains 
an open challenge, because dynamic measurements have 
to be resolved in both space and time 0] . 

Induced fluctuations are more easily accessible exper- 
imentally than spontaneous ones and can be related to 
one another by fluctuation-dissipation theorems. We in- 
troduce a dynamic susceptibility defined as the response 
of the correlator F(t) to a perturbing field x: 



Xx{t) 



dF(t) 
dx 



(1) 



The relaxation time of supercooled liquids increases 
abruptly upon cooling, so a relevant perturbing field 
is temperature, in which case Eq. ^ becomes xt(^) = 
dF(t)/dT. Density also plays a role in supercooled liq- 
uids, although a less crucial one Hence, another 
interesting susceptibility is Xp(t) — dF(t)/dP, where 
P is the pressure. Colloidal hard spheres undergo a 
glass transition [l9l | at high particle volume fraction 
tp. Thus, the appropriate susceptibility for colloids is 
X<p(f) = dF(t)/d<p. Equation n also applies in the fre- 
quency domain, Xx( w ) = dF(u>) / dx, where F{w) can be 
the dielectric susceptibility. We will show below that 
linear response formalism and fluctuation theory can be 
used to relate Xx(t) to the spontaneous fluctuations of 
C(i), and thus to Xi{t)- Thus, Xx(t) is an experimen- 
tally accessible multi-point dynamic susceptibility that 
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directly quantifies dynamic heterogeneity in glass form- 
ers. 

For molecular liquids, the dynamics conserves energy, 
volume and number N of particles, and one can estab- 
lish, in the NPT ensemble relevant for experiments, the 
following fluctuation-dissipation theorem, 

k B T 2 XT (t) = N(SC(t)SH(Q)) (2) 

where k B is the Boltzmann constant, H(t) the fluctu- 
ating enthalpy per particle, and C(t) the instantaneous 
value of a generic dynamic correlator F(t ). Both C(t) 
and H(t) are sums over local contributions [2(j, NC(t) = 
p Jd 3 fc(f,t) and NH(t) = p^Jk B c P T J d 3 fh(r,t). Here, 
p is the average number density, cp the constant pressure 
specific heat that sets the scale of the enthalpy fluctua- 
tions, (SH 2 ) = k B cpT 2 , so that the field h(f,t) has unit 
variance. Using translational invariance, Eq. |21 can be 
rewritten as: 

^T XT (t) = P J d 3 r(Sc(r, t)Sh(Q, 0)) (3) 

This expression shows that xt(^) directly probes the 
range of spatial correlations between local fluctuations 
of the dynamics and that of the enthalpy. In the case 
of colloids, the dynamics only conserves density and a 
similar expression can be obtained 

y/pk B TKT<px v (t) = P J d 3 r(Sc(r?,t)5p{d,0)) (4) 

where kt is the isothermal compressibility and 6p denotes 
density fluctuations rescaled by their root mean square. 

Equations |3 and 0| show that Xx(t) probes the extent 
of spatial dynamic correlations that differ from the ones 
studied in earlier theoretical and numerical works, which 
focused instead on x±{t) = 9 J d 3 r(Sc(r,t)Sc(0,t)). We 
have, however, established a direct relation between Xx{t) 
and Xi{t) by using the thermodynamic formalism de- 
veloped in [21|, which is generically applicable to bulk 
glass formers above the glass transition. For dynam- 
ics conserving energy and volume, we relate the fluc- 
tuations of C(t), and therefore Xi{t) measured in the 
NPT ensemble, to its isobaric-isoenthalpic counterpart, 
Xi v>)> which quantifies the amplitude of the fluctua- 
tions of C(t) in the NPH ensemble in which all configura- 
tions have exactly the same enthalpy: X &(t) = Xi FH (t) + 
k B T 2 x%(S)lcp. Because X4 PH (t) > 0, one derives an 
experimentally measurable rigorous lower bound | 22| for 
X4(i): 

X±{t) > —T 2 X 2 T {t) (5) 
cp 

A similar inequality holds between Xi{ u ) an d Xt(w), 
where X4( w ) denotes the amplitude of spontaneous fluc- 
tuations around F(u>). Similar arguments also apply to 
Xi{t) computed in the NVT ensemble preferred in nu- 
merical simulations. In that case, energy replaces en- 
thalpy in Eqs.[5|and|31 and the specific heat at constant 



volume, cy, replaces cp in Eq. [3] and relation|Sl Finally, 
we find that an inequality similar to relation [S] holds for 
colloidal systems, for which the volume and the number 
of particles are conserved quantities: 

X A{t) > P kBTKTw\%(t) (6) 

We have determined Xx(t) experimentally and numer- 
ically in three representative glass formers. For super- 
cooled glycerol near T g f=a 185 K, the real part of the 
dielectric susceptibility, e'(u>), was measured every 1 K 
in the temperature range from 192 K to 232 K. After 
fitting to a Havriliak-Negami form [[J, we use smoothed 
finite differences to evaluate xr(w) = d[e' (uj) / e' (0)] / dT 
and show in Fig. 1A the right hand side of relation 
as a function of inverse frequency. We plot in Fig. IB 
the right side of relation for hard sphere colloids where 
X<p(t) = 9f(q, t)/dip. The normalized intermediate scat- 
tering function |20j| f(q,t) is measured by dynamic light 
scattering [2j| for a wavevector q close to the first peak 
of the static structure factor. Several packing fractions 
are studied, from diluted samples where f(q, t) decays ex- 
ponentially in ~ 1 ms to concentrated suspensions with 
a two-step decay and a final relaxation time of ~ 10 s. 
Finite differences of data sets obtained for nearby ip are 
used to deduce XvW- Finally, we show in Fig. 1C nu- 
merical data obtained by standard molecular dynamics 
simulations of a binary Lennard-Jones mixture, a well- 
studied model for fragile supercooled liquids |24ll25l |. The 
dynamics is recorded at nearby temperatures through the 
self part of the intermediate scattering function, whose 
characteristic decay time s pan s a range from 1 ps to 100 
ns (using Argon units 

Dynamical susceptibilities behave similarly in all three 
cases. All display a peak for t w r Q , the average relax- 
ation time. The peak height increases when the glass 
transition is approached. This behavior represents the 
central result of our work. Together with Eqs. to 
it provides direct evidence of enhanced dynamic fluctua- 
tions and a growing dynamic lengthscalc associated with 
the glass transition. 

How tight the bounds of relations [S] and [S] arc depends 
upon the specific material and range of parameters stud- 
ied. A quantitative answer is given by simulations where 
the microcanonical quantity X4 llcr °(*)j i- e - the difference 
between Xi{t) m the NVT ensemble and kE.T 2 XT(t)/cv, 
can be easily measured. For the Lennard-Jones mix- 
ture, we find that the right side of relation [5] is much 
smaller than X4(0 a t high T, but the difference rapidly 
diminishes when T decreases. Both sides of relation [3] 
become comparable for the lowest temperature shown 
in Fig. 1C, which is still well above T g . Following 
Ref. |26l | , we also find that mode-coupling theory predicts 
Xi{t) ~ XtW ~ (T/T c — 1)~ 2 near the mode-coupling 
singularity T c > T gi provided that conserved variables 
are properly taken into account. 

These results support the idea that relation [S] can be 
used as an equality to quantitatively estimate X4^) a t 
low temperature, at least for fragile systems. This use 
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t (LJ units) 

FIG. 1: Dynamic susceptibilities in "x4 units", right side of 
relations |S| and |S| for three glass- formers. (A) Xt(oj) was ob- 
tained for 99.6% pure supercooled glycerol (in a desiccated 
Argon environment to prevent water absorption) by using 
standard capacitive dielectric measurements for 192 K < T < 
232 K (T g « 185 K). (B) XvtP) was obtained in colloidal hard 
spheres via dynamic light scattering. The static prefactor, 
pknTuT, was evaluated from the Carnahan-Starling equation 
of state |2(|. From left to right, tp = 0.18, 0.34, 0.42, 0.46, 
0.49, and 0.50. (C) xr{t) was obtained in a binary Lennard- 
Jones (LJ) mixture via numerical simulations. From left to 
right, T = 2.0, 1.0, 0.74, 0.6, 0.5 and 0.465 (in reduced LJ 
units lUI^). Relative errors at the peak are at most about 
10% for (A) and (C), and 30% for (B). For all of the systems, 
dynamic susceptibilities display a peak at the average relax- 
ation time whose height increases when the dynamics slows 
down, which is direct evidence of enhanced dynamic fluctua- 
tions and a growing dynamic length scale. 



of relation [3] is equivalent to assuming that dynamic 
heterogeneity in molecular liquids is strongly correlated 
with enthalpy fluctuations, and through a similar argu- 
ment, with density fluctuations in colloids. In fact, sup- 
posing that enthalpy is the only source of fluctuations, 
SC ~ (dC/dT)pSH/cp, and if we use the definition of 
cp, we obtain directly that Xi{t) ~ ^bT 2 x\^)I c p- A 
more general result can be obtained by taking into ac- 
count that energy and density are both fluctuating quan- 
tities, in which case X4 W is the sum of two contributions: 
Xi{t) « k B T 2 x 2 T {t)/cv + pk B Tn T p 2 x 2 p(t)- The second 
term is negligible in most fragile liquids |l8j . but dom- 
inates in colloidal systems. The presence of additional 
sources of fluctuations justifies that the rigorously de- 
rived inequality [S] docs not hold as an equality. 

Our results for dynamic fluctuations provide an esti- 
mate of the size £ of dynamic heterogeneity in liquids near 
T g . Because X4-W = pj d 3 r{5c(r,t)5c(0 7 t)), this quan- 
tity, once divided by the amplitude of the fluctuations 
at zero distance, (<5c 2 (0, t)}, defines a correlation vol- 
ume. The correlation functions are normalized to unity 
at t = 0, so (Sc 2 (0,t)) is of order one or smaller. Our 
simulations indeed show that in the temperature regime 
where the dynamics slows down and on timescales not 
much longer than the system relaxation time, this aver- 
age is of the order of one and displays extremely weak 
temperature dependence, as expected for a local quan- 
tity in glass formers. Thus, the height of the peak in the 
dynamic susceptibility, xt ~ ^b^IXt) 2 I°p-> yields di- 
rectly a correlation volume expressed in molecular units, 
x\ — (C/°)^j where a is the molecular size and cp is 
expressed in units oikp- Numerical 0,0 and theoreti- 
cal H [H [H HI] works suggest that < « 2 to 4. 

A direct comparison between our data and existing 
measurements can be performed for glycerol, where mul- 
tidimensional nuclear magnetic resonance (NMR ) exp er- 
iments show that £ = 1.3±0.5 nm for T = 199 K [27l li^ . 
Assuming a simple compact geometry for heterogeneities, 
C = 3, we estimate that £ increases from 0.9 nm at 
T = 232 K to 1.5 nm at 192 K. Given the assumptions 
involved in both approaches, and the uncertainty about 
numerical prefactors of order unity, the agreement is re- 
markable. An important physical conclusion of our work 
is that dynamic heterogeneity is strongly correlated to en- 
thalpy fluctuations in fragile liquids, although there is no 
signature of any static large-scale correlations 0, 0, 0] . 

For other glass- forming liquids, we obtain an estimate 
for £ at T g by assuming for simplicity that correlators 
obey time-temperature superposition, F(t) = J-(t/T a ), 
and using relational as an equality. One gets 

The logarithmic derivative is proportional to the well- 
known "steepness index" m, introduced in the glass lit- 
erature to characterize the fragility of glass-forming liq- 
uids |2^| . From reported values 0, |2^| of the quantities 
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FIG. 2: (A) Correlation volume xt(Tg) in supercooled liquids 
at the glass transition. Filled squares represent a lower bound 
to xt(Tg) in molecular units estimated through Eq. [7] Dif- 
ferent points represent different materials, which are ranked 
by their fragility m. Open squares represent the same quan- 
tity evaluated from available multidimensional NMR data ; us- 
ing |2|| xl = / d 3 rexp{-2r/CU forglycerol (T 9 + 10 K) j^, 
orthoterphenyl (T g + 9 K) pi l34| and d-sorbitol (T g + 7 
K) l28ll . A linear fit to the weak increase of xl with fragility 
is shown as a dashed line. (B) Correlation length £(T g ) in 
supercooled liquids at the glass transition expressed in bead 
units a. The correlation volume is first evaluated using Acp 
instead of cp in Eq.[7| Following l30ll . Acp is expressed in kn 
"per bead" units accounting for different molecular shapes 
and sizes. Using f = 3, the result is finally converted into 
a length scale expressed in bead units. The known empiri- 
cal correlations [2Sj between m, (3 and Acp translate into a 
weak increase of £(T 9 ) with fragility, which we fit with a linear 
relation shown as a dashed line. 



appearing in Eq. [7| and assuming a stretched exponen- 
tial form for T(x) = cxp(— x 13 ), we estimate x\{T g ) for 
different glass-forming liquids in Fig. 2A. 

For complex molecules, fluctuations that are unrelated 
to the glassy dynamics might contribute to the specific 
heat. These effects may be taken into account by re- 
placing cp in Eq. by Acp, the jump in specific heat 
at T g , which is sensitive only to the glassy degrees of 



freedom. Furthermore, for large molecules, the molecu- 
lar size is probably not the relevant microscopic length- 
scale, and it is sensible to express the specific heat in 
units of k a p er "independent bead" instead of molec- 
ular units |30| . These physical assumptions are used in 
Fig. 2B, where we have converted our results into length- 
scales expressed in bead units, and they lead to a trend 
similar to that of the main plot but with less scatter: Dy- 
namic correlations revealed by \t increase weakly with 
fragility [3l]|. This result is compatible with some theo- 
retical approach [3^J , but contrasts with others that pre- 
dict an opposite trend [3j| ■ This discrepancy might arise 
from the existence of at least two physically distinct dy- 
namic lcngthscales, one revealed by xt, and a second 
associated to xt- Although we found that both quan- 
tities are comparable for fragile systems, the bound in 
relation [3] may underestimate Xi f° r strong materials. 

To further test our length scale estimate of Eq. [7| we 
apply the formula to a polymeric liquid poly(vinyl ac- 
etate) (PVAc) using the monomer size for a [23 ■ We find 
£ « 2.0 nm at T g , to be compared with the value of 3.7±1 
nm obtained at T g +10K l^. |2§l [w e assume £ = 3, and use 
the data on PVAc given in |£}[2g]. Again, the agreement 
is satisfactory. A similar agreement is found for orthoter- 
phenyl and sorbitol, for which available NMR data are 
reported in Fig. 2A. Hence, we find that typical values for 
the dynamic correlation length at T g obtained via Eq. 
are in good agreement with previous experiments per- 



formed near T, 
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However, our approach 



has a broader scope, because it allows one to extend ex- 
perimental studies of dynamic heterogeneity to a range 
of temperatures not previously accessible and to the full 
time dependence of the fluctuations (Fig. 1). Finally, we 
remark that even for (strong) Arrhcnius molecular liquids 
with activation energy E, relation [3] and time tempera- 
ture superposition give x|(T) > (k B /c P ) x E 2 / (ksT) 2 , 
showing that dynamic heterogeneity must also exist in 
that case |35l |. in agreement with the general argument 
that for systems with finite range interactions, diverg- 
ing time scales must be accompanied by diverging length 
scales. 

Our experiments provide a quantitative demonstration 
that dynamic correlations and length scales increase as 
the glass transition is approached. More work is needed 
to characterize the time and temperature dependencies 
of dynamic fluctuations over a larger range of materials 
and parameters. Open issues also concern the precise 
space-time geometry of dynamic heterogeneity that fixes 
the value of the exponent £ and the relation between time 
scales and length scales, the connection between coopera- 
tivity and heterogeneity, and the extension of our results 
to the nonequilibrium aging dynamics encountered in the 
glass phase. 
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